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By introducing a boundary condition for the quantum wire, the Hubbard model is solved exactly 
by means of Bethe ansatz. The wave function for the bounded state is clearly defined, and the 
secular equation for the spectrum is exactly obtained. The ground state and low-lying excited 
states are studied in the thermodynamic limit. The ground-state energy in the strong coupling limit 
is obtained explicitly, and compared with the results of periodic boundary condition. 
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It is known that the behavior of one-dimensional elec- 
tronic system differs in many aspects from that of two- 
dimensional and three-dimensional systems M . The one- 
dimensional Hubbard model provides the opportunity to 
study correlation effects. There has been many discus- 
sions on Hubbard model since the Lieb-Wu solution was 
found This solution is an exact result with periodic 
boundary condition which can be thought of as on a ring. 
Owning to the dramatic achievement in nanotcchnology 
in recent years, the boundary effects of a sample will be- 
come more important. The Hubbard model on an open 
chain with completely confined ends was considered 
By adding some boundary field at the ends, the model 
was considered and discussed ||. The models |^-||] 
were easier and simpler as the wave function outside the 
wire is null. In present letter we study a more realistic 
model for quantum wires where the possibility of a non- 
vanishing wave function outside is taken into account. 
Using Bethe ansatz and some other methods in the lit- 
erature of integrable models, we obtain an exact result 
for the model hamiltonian. The ground state and the 
low-lying excitations are studied in the thermodynamic 
limit. The ground state energy in the strong coupling 
is shifted from the one for periodic boundary condition 
in a small amount. Both the spinon excitation and the 
holon-antiholon excitation are gapless. 

We consider a quantum wire of length L described by 
the hamiltonian 
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from either Hubbard ring or the Hubbard chain with 
reflection ends || or with boundary fields The Vj, 

and Vr in ([!]) represent the kinetic energies which an elec- 
tron must lose if it escapes out of the wire from left end or 
right end respectively. Instead of periodic or completely 
confined boundary condition, we will treat in this model 
with a more physical boundary condition. It is conve- 
nient to consider the states that span a Hilbert space of 
iV-particles 
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The eigenvalue problem H\if) >— E\if) > becomes an 
iV-particle quantum mechanical problem with the first 
quantized hamiltonian, 
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where Ajtp :— yb{- ■ ■ , xj + 1, • • •) + ip(- ■ ■ , Xj — 1, • • ■), 
V{xj) = V L 0(-Xj - L/2 - 1) + V R 6(xj - L/2 - 1), 6(x) 
the step function. The continuity limit of eq. (|^) concerns 
a equation in ref . . 

The number operator of electrons in the wire is de- 
fined as N — CiaCia, and that in the left side 

or the right side is given by N L = c ta C ™ or 

Nr — Cta,Cia respectively (summation also 

runs over spins a = 1/2,-1/2). The total number op- 
erator N t = N L + N R + N always gives N t \ip >= N\ip > 
in the Hilbert space of iV-particles. The boundary condi- 
tion which will be used to determine the spectrum reads 



lim C+C m |V->=0. 
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where Ci a is the operator annihilating an electron with 
spin component a on site i, and fii a := Cf a Ci a the local- 
number operator of electrons. Clearly, it is a Hubbard 
hamiltonian added with some terms which are defined 
outside of the wire. So the present model is different 



The Schrodinger operator (g) is invariant under any 
permutation of Sn, but is not invariant under transla- 
tion. Thus the total momentum of the system is not 
conserved, and the reflective waves must be taking into 
account. The wave function of Bethe ansatz form in the 



1 



region x G C(Q) := {x\ -L/2 < x Q1 < ■ ■ ■ < x QN < L/2} 
reads 

M*)= E MP, Q)e l{PklQx \ (4) 
Pew B 

where x := [xi,x%, ■ ■ ■ ,xn) with Xj G Z\ a := (aQi, 
a Q2, ■ ■ ■ ,0-qn), Q>j stands for the spin component of 
the jth particle; Pk (or Qx) is the image of a given 
k := (ki, &2, • • ■ , fcjv) (or x) by a mapping P G Wb (or 
Q G W A ); (Pk\Qx) = Ef=i( pk )j(Q^)j- The coeffi- 
cients A(P, Q) are functionals on Wb ® Wa- One may 
notice that the sum runs over the Weyl group 0] of the 
Lie algebra (denoted by Wb) but the wave function 
is defined on various Weyl chambers corresponding to the 
Weyl group of the Lie algebra An_i (denoted by Wa)- 

Any element of the Weyl group Wb can be ex- 
pressed as a product of the neighboring interchanges, 
a 3 : (• • ■ , zj, Zj+t, ■ ■ ■) i — ► (■ • • , Zj+i, Zj, ■ ■ •) and mirror 
reflection r 1 : [z\, Z2, '■•)*-* (—21, Z2, ■ ■ ■), where Zj de- 
notes either their image by a mapping. The 
requirement of anti-symmetry is (cr 3 ip) a (x) — —ip a (x), 
which gives 

A{P- a 3 Q) = -V Qj ' Q ^A(a j P; Q), 

where the spin labels are omitted and •p < 3-?»Q0'+ 1 ) is the 
spinor representation of the permutation a 3 . The Scat- 
tering matrix (S-matrix) which relates the coefficients 
A's between distinct regions in the configuration space 
of N electrons reads 

SQJ.QU+1) = sin(Pfc) J -sin(Pfc) 3+1 +zcPQ^+ 1 ) 
sin(P/s)j — sm(Pk)j+i + ic 

where 2c = U/t. Eq.(|) is the same as that in the 
Lieb-Wu solution [0. The coefficients A's in any re- 
gion are determined up to an overall factor by the 
§Qj,QU+V .- _7?Qj,QU+i)£<W,Qa+i) and the R L (or 
R R ), i.e., A(a 3 P;Q) = gQhQU+V A(P; Q), A(t 1 P-Q) = 
R L [(Pk)i]A(P;Q), A(t n P;Q) = R R [{Pk) N ]A{P-Q). 
The reflection matrices Rl and Rr are solved from the 
Schrodinger equation near the ends of the wire by taking 
into account the boundary condition ([}]), consequently, 

-ikL-i$ L (k) 
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R L {k) = 
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where iO(k) = \n[(e K+ik - 1) / (e K ~ lk - 1)} with 2t (cosh k- 
cos k) = Vl or Vr in 6l or Or . 

One may notice that the condition (||) does not give 
any constraints directly on ip but gives that on ip L and 
ip R . In deriving (^J) we have considered the wave func- 
tions outside of the wire, which actually vanishes in ||] 
but does not in present case. The wave function in the 
region xq\ < —L/2 < xq2 < ■ ■ ■ < xqn < L/2 takes 
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where Q'x = (0, XQ2, ■ • • , xqn) and Kj > with 
2t(coshnj — coskj) = Vl; in the region —L/2 < xqi < 

XQ2 < ■ ■ < %Q(N—1) < L/2 < x QN , 
N 

^R( X ) = J2 E A R {a"q'^Q)e- fi ^ N e l ^"^' k \ Q " x) 

j=l cr"eW" 

here Q"x = {xqi, ■ ■ • ,xqin~ 1), 0) and Kj > with 
2t(cosh Kj — coskj) = Vr. In the above, some nota- 
tion conventions are adopted , i.e., two subgroups of 
the Weyl group of B N , W' := {a 2 , a 3 , ■ ■ ■ , ct^" 1 , t n } 
and W" := {a\ a 2 , • • • , a N ~ 2 , a N ~ W" 1 }; 



two particular cycles <?'■' 



(xi, • • • ,Xj, ■ ■ ■ ,x N ) 



i-> (xi,- ■ ■ ,Xj-i,Xj + i,- ■ ■ ,x N ,Xj) and q'^ 

(xi,---,Xj,--- ,X N ) t-> (xj,Xl,--- ,Xj-i,Xj+i,- ■ ■ , X N ). 

As there are several identities in the Weyl groups, 
several consistency relations must be checked for the 
above solution. The Yang-Baxter equations || arising 
from both A(a 3 a 3+1 a 3 P;Q) = A(a 3+1 a 3 a 3+1 P;Q) and 
A(P;a 3 a 3+1 a 3 Q) = A(P;a 3+1 a 3 a 3+1 Q) are fulfilled 
identically. The Yang-Baxter equations with reflections 
arising from A{T 1 a 1 T 1 <j 1 P;Q) = A(a 1 T 1 a 1 T 1 P;Q) and 
A(T N a N - 1 T N a N ~ 1 P;Q) = A((t n ~ 1 t n a N ~ 1 r N P; Q) 
are satisfied identically because the reflection matrices 
are just scalar factors. 

As either {a 1 , ■ ■ ■ , o^" 1 , r 1 } or {a 1 , ■ ■ ■ , a N '\ t n } 
can be chosen as the basic elements of Weyl group Wb , 
we must consider another consistency condition for Rl 
and R R . Using r 1 = uV 2 • • • a N - 1 T N a N ~ 1 ■ ■ ■ a 2 a x and 
applying the S'-matrices successively, we obtain an eigen- 
value equation in the spinor space: S 1,2 {r\\ +r)2)S 2,3 (r]i + 

m )S 1 ' 2 ( Vl - m)A(P;Q) - RLKPk^R^KPk)!] 
A(P; Q), which is equivalent to 



S^im + ^)5 1 ^(? 7l +%)••• S^( m + r) N ) 
S 1 ' N (m-m)---S 1 ^(r, 1 -r]2)A(P;Q) 
= R L [(Pk) 1 ]R B 1 [(Pk) 1 ]A(P;Q), 
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where r\j = sin(Pfc) 3 . This relation guarantees 
the consistency for any reflection (ki, ■ ■ ■ , kj, ■ ■ ■) — > 
(fci, • • • , — kj, ■ ■ ■)■ 

The eigenvalue problem (Q) can be diagonalized by 
means of Sklyanin approach H, in which the transfer 
matrix is defined as t(a) = trR' L 1 (a)T + (a)RR(a)T^(a) 
where (a) = T AN (aTaN) • • ■ T AN (aTa2)TAN(a : Foti), 
and T Ai (a) := S" Aj '(a) G End(V A ®V 3 ). The Bethe 
ansatz equations are obtained as follows 
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where Sg(x) = (x — i(ic)/(x + i(3c). A set of coupled 
transcendental equations are derived by taking the loga- 
rithm, 

*i = (£)4-^[0£(*i)+M**)] 

M 



[©i(A„ - A„) + 9i(A^ + A,)] = -2ttJ, 

i/=i 
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where 0/3 (x) := 2tan (x//3c), — 7r < tan _1 (x) < 7r. The 
quantum numbers of both the quasi-charge and quasi- 
spin, Ij and J M , take integer values regardless of N — M 
being even or odd. Eq.(||) becomes the result of Ref. ||] 
in the limit Vl,r — ► oo. 

We consider the thermodynamic limit by introducing 
the densities (distribution function) of roots, p(k) and 
ct(A) , and that of holes Ph{k) and 07i(A). Then the sec- 
ular equations (|j) become the following coupled integral 
equations 



p(k)+ Ph (k) 
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+ cosfc / dA'Ki(smk\\')a(\'), 
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where Kp(x) := 7r~ 1 /3c/(/3 2 c 2 + x 2 ) and a notation 
F(x\y) := [F(x — y) + F(x + y)]/2 is used here and from 
now on. The B and D are determined from the condi- 
tions 
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where the factor 1/2 is necessary [[To| 

The ground state of the present model is a Fermi sea 
described by po(k) and cto(A), where po(k) is the distri- 
bution function of charge with momentum k and (To (A) 
that of down spins with respect to the rapidity A. The 
distributions of the roots satisfy ([)]) with p h — 6(0) /L, 
<Jh = and B = 00. The energy of the ground state can 
be calculated once po(k) is known. In this case, the dis- 
tribution functions are solved in a closed form by Fourier 
transform, 
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where we have used the following definitions, 
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Obviously, S\ (x) is the conventional hyperbolic function 
sech(a;) which occurs in Lieb-Wu solution. 

Eqs.(|i"l|) can be solved explicitly in some special cases. 
In the strong coupling limit c 3> 1 , it is simplified to 

2n 2n An N 

27rp 8tr (k) = 2+ —p h (k) - —6(0) + — (— )ln2cosfc, 
L L c L 

(12) 



where 2-Kp\ ) (k) = — + 9n(k)) arises from the 

dk 

boundary contribution at the ends of the wire. The Do is 
completely determined from the explicit form ( |T2] ) , hence 
the energy density of the ground state is obtained 
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where 



£fc = 
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Tj(fc)] sin /cdfc 
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which vanishes when Vl.r. 00 ■ It is worthwhile to com- 
pare (|l^) with the result in the case of periodic boundary 
condition. The difference is calculated, 



E - E p e 
L 



where Eq CT stands for the ground-state energy with pe- 
riodic boundary condition. Clearly, the difference is rel- 
atively small for a large system which agrees with the 
common understanding about the effects of boundary 
condition. 

It is convenient to study the excitations by introducing 
p(k) = p (k)+pi(k)/L andcr(A) = ao(\)+ai(\)/L where 
Po(k) and cto(A) satisfying (H) with Ph(k) = 6(0) /L and 
Ch(A) = 0. The excited energy up to the order 0(1/ L) is 
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where p stands for the chemical potential |TT[. The jl^ ) 
is related to pi(fc) only and is valid for both the spin 
excitation and the charge excitation. For a large system 
the D can be replaced by Do- 

The simplest spin excitation is a triplet, i.e., two-hole 
state with 07, (A) = [6(X | A x ) + S(X \ X 2 ))/L. The 
excitation energy is composed of two terms AE tr i P = 
£s(Ai)+£ s (A2), each of them can be identified as a spinon 
excitation energy 
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where pf(k, X) solves the following integral equation, 
cos k 



pl(k;X) 
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From the asymptotic behavior of the Si one can find that 
the spinon excitation is gapless. 

The excitation in charge configuration is a variation 
from the ground state in the sequence of charge quantum 
numbers {Ij}. The simplest case is the holon-antiholon 
excitation, i.e., one hole kh € [— Dq,Dq] and one 'parti- 
cle' k p [— D ,D ] state. In this case pt = S(k | kh)/L, 
the integral equations @ give 

/>oc 

pl(k) = -S(k | k h )+cosk / dX'Ki(smk\X')a1(X'), 



a c i(X) = Ki{X\smk p ) + / dk'Ki (A| sin k')p c x {k') 
J -D 



dX'Ki{X\X')(7t{X'). 
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The excitation energy is composed of two terms AE = 
£c(kh) ~ £c(k p ). They can be identified as a holon ex- 
citation energy s c (kh) and an antiholon (particle state) 
excitation energy e c = —e c {k p ) respectively, namely, 



e c (k) = 2t cos k 



dk(2tcosk +fi)pl(k,k), (18) 
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where the p\{k 1 k) is determined by 
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Since both kh and k p can tend to Dq, the holon-antiholon 
excitation is of no gap. Further case is the holon-holon 
excitation which involves complex fc-pairs and is gapful. 

The excitation in spin configuration is allowed to form 
a singlet, i.e., two holes and a 2-string A ± = (Ai+A2)/2± 
ic/2. The excitation energy is composed of two terms 



AE sing = e s (Ai) + e s (X2). 



(19) 



These two terms are free spinon excitation energies given 
by (HHH). 

In summary, we obtained an exact solution for the 
Hubbard model with the realistic quantum-wire bound- 
ary condition that the possibility of a nonvanishing wave 
function outside is taken into account. The ground state 
and some low-lying excitations are studied in the ther- 
modynamic limit. 
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